Splitting of roton minimum in the z/ = 5/2 Moore- Read state 
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We calculate the dynamical structure factor of the u = 5/2 non-abelian quantum Hall state in 
the dipole approximation, valid for large momenta. Due to the fact that both quasi-particles (qps) 
and quasi-holes (qhs) have an internal Majorana degree of freedom, a qp-qh pair has a fermionic 
degree of freedom which can be either empty or occupied, and leads to a splitting of the roton 
mode. Observation of this splitting by means of finite wavelength optical spectroscopy could provide 
evidence for Majorana modes in the p — 5/2 quantum Elall state. 
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The J/ = 5/2 quantum Hall (QH) state is expected to 
support non-abelian quasi-particles (qps), which could 
be used for topological quantum computation There 
is both theoretical and experimental evidence that the 
1/ = 5/2 QH state is indeed described by the non-abelian 
Moore- Read state Q . Most significantly, the quasiparti- 
cles in the Moore- Read state have charge e/4. This has 
been confirmed by tunneling shot noise Q and in- 
terference [Hj experiments. However, photoluminescence 
experiments seem to indicate that the state is not spin 
polarized, contrary to the Moore- Read state @. Recent 
resonant light scattering experiments suggest that both 
polarized and unpolarized domains may form in the sec- 
ond Landau level in general Q . Clearly more evidence is 
needed to confirm the validity of the Moore-Read state. 

The roton mode is the collective excitation of quantum 
Hall states In integer quantum Hall systems, it can 

be understood as a particle hole pair, in fractional quan- 
tum Hall (FQH) systems as a quasi-particle (qp) quasi- 
holc (qh) pair. For hierarchical FQH states, there is one 
roton minimum for each level of the hierarchy. The roton 
mode can be probed with the help of optical excitations, 
but only in the presence of finite wavelength density mod- 
ulations. Using surface acoustic waves to create a charge 
density modulation, several roton minima for hierarchical 
FQH states were recently observed experimentally [lo|. 

The non-abelian statistics of = 5/2 qps is encoded by 
an internal Majorana degree of freedom. Here we show 
that as a direct consequence of these Majorana fermions, 
the roton mode in the i/ = 5/2 Moore- Read state splits 
into two branches. The two Majorana modes of a qp-qh 
pair can be combined into a neutral fermion, which can 
can be either occupied or unoccupied. Due to the finite 
overlap of the two Majorana wave functions associated 
with qp and qh, an occupied neutral fermion has higher 
energy than an unoccupied one, hence the splitting of 
the roton mode. The splitting both oscillates and decays 
exponentially with increasing separation of the pair. The 
main result of this work is the excitation spectrum given 
by Eq. (fTO|). whose form is shown in Fig. 1. We argue that 
the splitting of the roton mode and hence the non-abelian 




FIG. 1: (Color online) The expected roton dispersion. The 
thick lines are given by Eq. (|10p . constituting the dipole ap- 
proximation, valid f or q beyond the roton minimum, which oc- 
curs near qlo ~ 1.4 [IJl. The usual single roton mode (shown 
with black dots) has split into two branches with distinct ex- 
trema, denoted lo+ and oj-, where a;+ — lj- = 2E{2lQq). u- 
is the brighter primary branch, and is the secondary split- 
off branch. In the the small g-regime, the dipole Hamilto- 
nian is no longer valid and the roton minimum is taken to be 
parabolic with appropriate parameters used as outlined in the 
text. This regime is plotted with thin lines. A ~ 0.025e^/eio 
is the u = 5/2 exciton gap 15]. The grey shaded region de- 
notes the two-roton continuum. 



nature of the = 5/2 state can be directly probed in an 
optical experiment similar to [lo| . 

We model the Moore-Read state as a p-wave super- 
conductor of composite fermions and then use a 
Ginzburg-Landau approach to describe the dynamics of 
the superconducting order parameter. Charge e/4 qps 
of the 5/2 state are vortices of the superconductor, and 
are treated in the same way as the vortices of condensed 
composite bosons [13, in order to calculate the ro- 
ton branch. The dynamics of the order parameter effec- 
tive theory with N vortices are identical to those of an 
N point particle Hamiltonian with pairwise interactions 
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and a constant kinetic term per particle, with the com- 
mutation relations [xi,yj] = —it^Si^j/v, where i denotes 
the i*^ vortex, and /q = yJh/eB is the magnetic length. 

Vortex configurations of the order parameter in the 
Bogolubov-de-Gennes (BdG) Hamiltonian give rise to 
low-energy Majorana excitations, pairs of which form 
neutral fermions. Our Hilbert space then includes both 
the vortex degrees of freedom as well as the occupancy 
of neutral fermions that reside in pairs of vortices. We 
shall truncate the many-body Hilbert space of vortices to 
a single vortex-antivortex pair, precisely the qp-qh pair 
that forms the roton mode. 

For Majoranas separated by r, the unoccupied and oc- 
cupied energies of the neutral fermion have been calcu- 
lated for a p-wave superconductor [3| 



E{r) = ± 



2Ao cosippr + 7r/4) 



-exp 



(1) 



in which Aq is the mean field value of the order parame- 
ter field, pF is the Fermi velocity, and ^ is the Majorana 
coherence length. To make quantitative predictions for 



the Moore- Read state, we use the estimates of Ref. [17 1 
for the magnitude of the splitting near the roton mini- 
mum (w O.Ole^/e/o), the coherence length, f « 2.3Zo, and 
for pp « 7r/3/o- 

The Hamiltonian for a single qp-qh pair in this effective 
picture is 



Hair) = A + V{f) + E{f){2c^c - 1) 



(2) 



Where A is the energy gap to create two isolated qps 
[l5| . r = |ri — is the dipole separation operator, 
V is the Coulomb interaction between the charge ±e/4 
vortices, and = 71 — 172 is the neutral fermion cre- 
ation operator which contains two Majorana operators 
(71 and 72). The non- vanishing commutators in relative 
separation (f) and centre of mass (R) coordinates are 
[fx.Ry] = 2ilQ, and [fyjRx] = — 2iZg. The canonical cen- 
tre of mass momentum is then related to the relative sep- 
aration as f = 21qZ X K, where z is the out-of-pane unit 
vector. The Hamiltonian has the dynamics of a dipole 
drifting with constant velocity in a magnetic field, with 
an internal 2-level degree of freedom. 

For qp-qh separations comparable to the intrinsic size 
of a vortex, the dipole approximation breaks down, and 
the leading contribution to the dynamic structure factor 
is of quadrupole type. For separations much larger than 
the intrinsic size of a vortex, the single dipole approxi- 
mation is asymptotically exact [l3| . 

Formally, the ground state of the Hamiltonian Eq. ([2]) 
is at r = 0, with both a divergent splitting energy Eq. ^ 
and a divergent potential V. These divergences are an 
artefact of the point-particle description of vortices: the 
Hamiltonian Eq. ^ is only valid for well separated qps, 
ie. r greater than twice the intrinsic size of a single vortex. 
However, in the following calculation we use a completely 



overlapping qp-qh pair as the vacuum state and assign 
the energy zero to it. Then, the density operator de- 
fined below separates the qp and qh from each other and 
creates a finite energy excitation described by Eq. ([2]). 

The parity of the number of quasiparticles must be 
conserved: an excitation breaks a Cooper pair, necessar- 
ily incorporating two fermions. It is necessary then, to 
add to the Hamiltonian Eq. ([2]) a set of vortices pinned by 
disorder, which can, in pairs, support neutral fermions. 
These vortices will inevitably be present in any real sys- 
tem due to disorder, but their density can also be tuned 
externally. 

The arrangement of pinned vortices is assumed to be 
random. We are interested in the dynamics of a single qp- 
qh pair, created by a momentum fluctuation within this 
background of pinned vortices. If the average separation 
of two pinned vortices is greater than a few times the qp- 
qh separation, the overlap of the dipole vortices with the 
disorder vortices will be appreciable for at most one pair 
of disorder sites. Thus, in the spirit of the random singlet 
phase [l8|, we can discard all other pinned vortices, and 
reduce our system to that of four vortices: the qp-qh pair 
and the two closest pinned vortices. Let the operators c 
and correspond to the dipole neutral fermion, and b 
and to the single pinned vortex pair. We can now 
introduce a tunneling term to the Hamiltonian, 



tdbc ~ e^c^b^ + t'^Pc - t'J:c'ib + {2Pb - l)e<j (3) 



Here, ta and are the complex overlap integrals between 
the four sites (1,2,3,4), and denotes the energy split- 
ting of the localized neutral fermion. We obtain td = 
£;(|ri-r3|)-£;(|r2-r4|) + i(£;(|r2-r3|) + i;(|ri-r4|)), 
where t'^ has a similar form but is irrelevant here. We can 
now proceed to calculate the collective excitation spec- 
trum, given by = / uiS{q,u!)du}/ J S{q,u})duj 0], 
where the dynamic structure factor is given by 



S{q,uj) 



dte-*(p,(i)/5-,(0)). 



(4) 



Here Pq{t) is the charge density operator with time de- 
pendence /9g(i) = e^^* pq{0)e~^^* , in which = Hq+Ht- 
For a single qp-qh pair, Pq{0) = J2s=±i se*i (^''+^'. We 
calculate the eigenstates of Hq explicitly, and perturba- 
tively introduce Ht up to second order. Fourier trans- 
forming, we obtain two terms, S'o(g,a;) which contains 
only the unperturbed Hamiltonian, and Si{q,uj) which 
contains Ht to second order. The first term produces the 
primary excitation branch in which the neutral fermions 
are unoccupied. That the neutral fermions are unoccu- 
pied is immediately obvious as the charge density oper- 
ator acts only on the vortex sector of the ground state 
Hamiltonian. The second term contains only those terms 
cx tdt^, as the mixed creation-annihilation operator terms 
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vanish. This term in the structure factor contains opera- 
tors that fill the two neutral fermion states, and produces 
the secondary split-off branch, peculiar to the v — "0/2 
Moore-Read state. 

The total Hilbert space is a product of the vortex 
space with matrix elements of the form (0|/5g|K), and 
the neutral fermion space with matrix elements of the 
form Jq dti Jg'^ dt2{0\HT(ti)HT{t2)\m,n), where m and 
n correspond to the occupancy of the disorder-pair and 
the dipole respectively, and the time dependence of Ht {t) 
arises from using the interaction picture with respect to 
Hq. The unusual temporal limits reflect that for times 
^1(2) < or ti(^2) > the system is in the vacuum state 
and so cannot support Majoranas in the dipole, making 
Ht nonzero only in the time interval [0,t]. The vortex 
part of the ground state Hamiltonian is diagonal in rel- 
ative separation (f), and the centre of mass coordinates 
do not appear. The wavefunction can be written as 



iK R 



2l^z X K) 



(5) 



This is a point particle description of the vortices. How- 
ever, in order to calculate the overlap of two vortices, the 
vortices must be given a spatial profile whose width re- 
flects the size of the vortex. We introduce a Gaussian pro- 
file, such that S{r + 2l^zxK) 2?o/v^e-l'"+2'o«xK| V4/^ . 
For a charge-density fluctuation of momentum qx, we 
obtain a single solution which is a qp-qh pair with sep- 
aration 2lQqy, and matrix element 2 exp(— thus 
obtaining 



m.Ti 



(6) 



•|(0,0,0|pq|m,n,K)|2 

where the eigenvalues are 

e„„K - ^ + V{2llK) + {2n~l)E{2llK) + {2m-l)ed (7) 

and td is the contribution from the pinned vortex Majo- 
ranas. As the density operator is simply the real space 
translation operator, the summation in the Majorana sec- 
tor gives (5m,o^n,o- Therefore this part of the structure 
factor is the primary branch of Eq. without the de- 
creased spectral weighting of Mq. As for Si{q,t), the 
first order expansion of Ht gives zero, and for the sec- 
ond order term we find 



-i(em7iK — eooK)(il— 12) 



-i(eooK-eooo-l-<^)t|^0,0,0|pgC6|TO,n,K)|^ 



(8) 



Here the summation over the Majorana sector results in a 
single occupied neutral fermion term ^m,i^n,i due to the 



operator pair cb and its complex conjugate. Therefore 
this term contributes the delta functions in Eq. ([9|). We 
thus obtain the structure factor, and therefore excitation 
spectrum, of our qp-qh pair system where the vortices 
support neutral fermions with non-zero overlap integrals, 
on top of a sparse background of qps and qhs. 

The final form of the structure factor up to second 
order in Ht is 



S{Q,Lj) 



{1 - All)d{Lo^{Q) - uj) 



+ M^{6{uj+{Q) E{Q)6'{uj^{Q) - uj)) 



(9) 



Where Q = 2Zgg, and Mq = \td\/AE{Q), and 5' denotes 
the derivative of the delta function with respect to its 
total argument. This latter term effectively doubles the 
matrix element of the secondary branch, such that the 
spectral weighting of the higher energy branch is 2Mq, 
and the lower branch has spectral weighting 1 — 2Mq. 
Solving for the excitation spectrum, we obtain two dis- 
tinct branches given by 



uj±{Q) = A + V{Q)±E{Q) 



(10) 



Here, the splitting of the localized fermion has been 
neglected because in the relevant parameter regime it 
will be much smaller than the E{Q) due to the relatively 
large distance between localized vortices. The splitting 
of the two excitation branches is a direct result of vor- 
tices supporting Majoranas, and the exponentially decay- 
ing splitting demonstrates that well separated Majoranas 
have approximately zero energy. These are two necessary 
ingredients for the vortices to be non-abelian quasiparti- 
cles. The expected dispersion is shown in Fig. 1. Exper- 
imentally, the most important signature is the crossing 
of the two roton branches near qig = 2. As a conse- 
quence, the roton minimum has split into two and shifted 
to higher and lower q values, respectively. The param- 
eters used for E{Q) are outlined below Eq. ([T|). We 
have also taken the roton gap to be half the exciton gap 
A « 0.025e^/eZo, the energy required to create a maxi- 
mally separated qp-qh pair |19|. which agrees well with 
exact diagonalization studies [20|, |21|, and have set the 
roton minimum to occur at qlo ~ 1.4, in good agreement 
with numerics [l3 |. 

The two quantities E{Q) and td are independently tun- 
able. The former determines the amount of splitting be- 
tween the two branches and is a function of the separation 
of the dipole (which is proportional to the momentum at 
which the roton branch is probed). The latter, td, de- 
termines the spectral weighting of the secondary branch, 
but will also broaden both branches, and is a function 
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of the proximity of the dipole to the closest pinned vor- 
tex pair, and therefore is directly related to the disorder 
density. 

The parameter td must be tuned such that td <C E{Q), 
so that the broadening of the excitation spectrum lines 
does not fuse the two branches. However, td must be 
sufficiently large that a measurably bright higher energy 
band, whose intensity is proportional to \td\'^ is observed. 

In addition to the Moore-Read state, alternative wave- 
functions have been proposed for the description of the 
1^ — 5/2 state. Most notably, the 331 candidate wave- 
function differs from the Moore-Read state in that the 
spins are unpolarized. We can obtain information about 
the roton excitation of the 331-state by using the fact that 
this wavefunction has been used to describe pseudo-spin 
correlations in Z/'tot — 1/2 bilayer QH systems. For bi- 
layer systems, the inter-layer Coulomb interaction takes 
the form V(r) oc (r^ -I- d^)^/^, where d is the interlayer 
spacing. The relevant case for h' — 5/2 then, is to take 
the (i — >■ limit of these results. Two magnetoroton spec- 
tra have been predicted for i^tot = 1/2, due to in-phase 
and out-of-phase pseudo-spin density fluctuations [j^ . 
and these persist in the d limit. If the 331 wave- 
function correctly describes v — 5/2, then the following 
observations can be made. The splitting between the two 
branches is a minimum at the roton minima, contrary to 
our result. Depending on the splitting magnitude, the 
upper branch may be lost in the excitation continuum 
altogether. Furthermore, the splitting will not decrease 
exponentially (or even at all) with increasing q. Finally, 
the splitting will not oscillate. The combination of these 
signatures allows for a clear distinction from the split- 
ting of the roton minimum for the Moore-Read state as 
described above. 

Secondly, a 'generalized composite fermion' picture of 
1^ — 5/2 was recently discussed and illustrated by exact 
diagonalization In order to understand the lowest 

energy excitations, the possibility of two types of ele- 
mentary excitations was assumed: a magnetoroton con- 
sisting of a qp-qh pair, and a broken composite fermion 
pair. This scenario for v — 5/2 was raised already some 
time ago [2l|. From the numerical work [g^l, it seems 
possible that the gap for a broken pair type excitation 
may be slightly below the magnetoroton gap. However, 
the broken pair spectrum is a Bogolubov spectrum that 
monotonically increases with increasing q, while the mag- 
netoroton spectrum on the other hand asymptotically ap- 
proaches the exciton gap with increasing momentum. For 
momenta relevant to our current calculation then, that is 
for q beyond the roton minimum, the broken pair energy 
will be larger than the magnetoroton energy, and will 
disappear into the continuum region. Most importantly, 
composite fermions are charge neutral and are expected 
to couple only weakly to the charge density operator. For 
this reason, it is unclear whether the broken pair mode 
can be detected experimentally when studying the struc- 



ture factor S{q,U!). 

Finally, we comment on the neutral fermion gap: 
the energy required to add a single unpaired composite 
fermion to the system [2^ 24 1, which is comparable to 



the exciton gap A |23|. The neutral fermion dispersion 
also contains roton-like minima ^2^], however an exper- 
imental probe must change the particle number parity 
24{ . Thus the collective excitation spectrum presented 
here will not include the neutral fermion dispersion. 

We have shown that the magnetoroton dispersion of 
the z/ = 5/2 Moore-Read state splits into two distinct 
branches due to a fermionic degree associated with qp- 
qh pairs. This neutral fermion consists of a pair of Ma- 
jorana states in the vortex cores which are responsible 
for the non-Abelian statistics of the Moore-Read state. 
In order to overcome the requirement of parity conser- 
vation, we coupled the qp-qh pairs to localized qps in 
the bulk of the quantum Hall state. The intensity of 
the secondary branch of the roton dispersion is propor- 
tional to the square of the coupling strength to these bulk 
qps. Experimental observation of this splitting by means 
of optical spectroscopy would constitute compelling evi- 
dence for the validity of the Moore-Read state and thus 
non-Abelian quasiparticles in the v = 5/2 QH state. 
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